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Abstract 

We analyze the role of nuclear modifications of parton distributions, notably, the nu- 
clear shadowing and antishadowing corrections, in production of lepton pairs from decays 
of neutral Z and 7* gauge bosons in proton-lead and lead-collisions at the LHC. Using the 
Collins-Soper-Sterman resummation formalism that we extended to the case of nuclear par- 
ton distributions, we observed a direct correlation between the predicted behavior of the 
transverse momentum and rapidity distributions of the produced vector bosons and the pat- 
tern of quark and gluon nuclear modifications. This makes production of Zj^j* in pA and AA 
collisions at the LHC a useful tool for constraining nuclear PDFs in the small- a; shadowing 
and moderate- a; antishadowing regions. 
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1 Introduction 



Nuclear modifications of structure functions and parton distribution functions (PDFs) are a 
firmly established phenomenon extensively studied during last three decades [U [2J O [H E] . 
Presenting the emerging pattern of these modifications in terms of the ratio of the nuclear 
to free nucleon structure functions, R(F 2 ) = F2a(x)/[AF 2 n (x)], one observes the following 
characteristic trend in different regions of Bjorken x: R(F2) < 1 for x < 0.05 — 0.1 (nuclear 
shadowing), R(F 2 ) > 1 for 0.1 < x < 0.3 (antishadowing) , R(F 2 ) < 1 for 0.3 < x < 0.8 (the 
EMC effect), and R(F 2 ) > 1 for x > 0.8 (Fermi motion). The four nuclear effects that are 
responsible for the medium modifications of the PDFs have different magnitudes and patterns. 
Indeed, the suppression due to nuclear shadowing increases with the atomic number A and could 
be as large as, e.g., 20% for 40 Ca; the enhancement due to antishadowing is a few-percent effect 
that does not reveal a distinctive A dependence; the EMC effect suppression increases with A 
and reaches, e.g., 20% around x = 0.65 for 56 Fe; the enhancement of R(F 2 ) due to the nucleon 
Fermi motion inside nuclei for x > 0.8 increases with x, becomes as large as 20 — 30% for x ~ 0.9, 
and leads to a nonvanishing nuclear PDF at x > 1. 

In the leading- twist approach, this pattern of nuclear modifications of R(F 2 ) translates into 
a similar picture of medium modifications of PDFs [71 [5"1 [9| ITUl 111]. However, the magnitude of 
nuclear modifications of individual quark and gluon distributions in nuclei in different regions 
of the light-cone fraction x is rather uncertain. This is especially true for the gluon channel, as 
the indirect extraction of the gluon distribution from mostly fixed-target experiments on deep 
inelastic scattering (DIS) with nuclear targets leaves the gluon distribution in nuclei largely 
unconstrained. 

In this paper we examine the possibility of establishing additional constraints on nuclear 
PDFs, primarily in the small- x nuclear shadowing region, by analyzing the size of nuclear cor- 
rections in neutral gauge boson production in proton-lead and lead-lead reactions at the Large 
Hadron Collider (LHC). In particular, we study the transverse momentum distribution of Z/-f* 
production with its subsequent decay into a lepton pair in the reaction A+B — > (Z/j* — > ££)+X, 
with A = p or 208 Pb, B = 208 Pb, and i = e or fi, by using the Collins-Soper-Sterman (CSS) 
formalism [12 1 113 1 PHI I15| implemented in the resummation computer code ResBos |16lll7j. For 
this computation, we used the nuclear PDFs described in Sec. [3] and carried out the resummation 
calculation with the resulting nuclear parton distributions. The choice of this channel presents 
several advantages: it is very clean, has a high rate and is very well understood in proton-proton 
collisions. Moreover, the Qt distribution of heavy vector bosons can be easily reconstructed 
experimentally and, thus, can be measured very precisely. Access to a wide range of boson 
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virtualities provides a powerful test of the leading-twist factorization and QCD evolution. 

The extension of the CSS resummation formalism to nuclear scattering is described in Sec. [2j 
The treatment of nuclear corrections in the leading-twist approach is reviewed in Sec. El while 
Sec. U] presents numerical predictions for nuclear corrections in Z/~f* production in pA and AA 
scattering at the LHC, in two distinct intervals of the Drell-Yan pair mass, 70 < Q < 110 
GeV and 5 < Q < 20 GeV. We examine nuclear corrections to the Drell-Yan rapidity and 
transverse momentum distributions and discuss dependence on the nonperturbative Ht smearing 
corrections. Finally, Sec. [5] contains our conclusions. 

2 Application of Qt resummation to nuclear scattering 
2.1 pp collisions 

In the standard framework |12t [T3| [T^l [15] for transverse momentum resummation in nucleon- 
nucleon collisions, the differential cross section for vector boson production reads 



where Q, Qt, and y are the invariant mass, transverse momentum, and rapidity of the vector 
boson, respectively. The cross section consists of two terms, a Fourier-Bessel integral over 
the transverse position (6) that evaluates the resummed logarithmic contributions in the limit 
Q T <C Q 2 , and the finite piece Y(Q,Qx,y) that dominates at Q\ Q 2 . The resummed 
form factor W(b,Q,y) receives leading-power and power-suppressed contributions, and those 
are incorporated into the resummation integral of Eq. ([T|) as W pert and W NP , respectively. 

The form factor W pert and non-singular contribution Y (defined as the difference between 
the fixed-order calculation and the approximation for the perturbative part W pert to the same 
order in a s ) are computed in perturbative QCD. The explicit expression for W pert is 





Jo(q T b)W pert (K,Q, y)W NP (b, Q, y) + Y(Q, Q T , y) 



(1) 




j=u,d,s... 
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in terms of the Sudakov exponential 
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and convolutions \Cj/ a ® fa/h] of Wilson coefficient functions Cj/ a {z, C1/C2, Hf = C3/6) with 
distributions f a /h( z i I^f) of partons a inside the initial-state hadron /i. The normalization factor 
is 

(o) _ Mfc^ g 2 



12iV c 5 (Q 2 - Mf) 2 + M|r| 
1/2 



'l-4|e i |sin 2 6' W) ) +1 



sin 2 9 W ) + sin 4 



in — > [Z — > e + e )X, and 
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in pp — > (7* — > e + e~)X. Here is the Z boson mass, Gf is the Fermi constant, 9 W is the 
weak angle, ^/S is the collision energy, N c = 3 is the number of colors, chem is the fine structure 
constant, and e-,=2/3 or -1/3 is the fractional quark charge. The scale dependence of W v 
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ascribed to three scale parameters C\ = bfi, C2 = n/Q, where fi is the integration variable in 
the Sudakov exponent, and C3 = fipb, where \xf is the factorization scale in the convolutions 
[C j/a ® fa/h] with^the PDFs. 

The function W NP (b, Q) in Eq. (TjQ) provides a model for power-suppressed terms in the 
resummed form factor that play an important role at transverse positions of order 1 GeV -1 
or above. We evaluate the perturbative part W pert as a function of 6* = bj yl + b 2 /b^ax 
and combine with the nonperturbative function (b, Q). The 6* variable [13|. I14j realizes a 
smooth transition from the leading-power terms that dominate at small b to nonperturbative 
dynamics that is important at large b. We implement the 6* convention in its extended version 
[15] that can simulate various trends in the behavior of W(b, Q) at b > 1 GeV -1 suggested by 
the nonperturbative models [191 1201 Ell [221 [231 IS], while at the same time retaining the exact 
perturbative prediction for W(b,Q) at 6 < 1 GeV -1 . 

Parameters of W (b,Q) are constrained by a fit to the experimental data. To reproduce 
the observed Qt distributions in Drell-Yan process, it suffices to approximate W N (b,Q) by 
a Gaussian smearing factor exp [— a(Q, x%, X2)b 2 ~\ , where the b 2 dependence is associated with 
the lowest power-suppressed contribution [25J, and a(Q, x±, X2) is found from the fit [261 E]- A 
recent parametrization for a(Q) from [18] describes well both the low-Q pair and Z production. 
Recently, constraints on a{Q) at Q w Mz were updated [271 [28] using data on angular distribu- 
tions da/d(f>* in Z boson production at the Tevatron Run-2 [29]. This updated nonperturbative 
parametrization will be used in our simulations, as it provides excellent description of Z boson 
Qt distributions in pp collisions at the LHC. 
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2.2 pA and A A collisions 

The CSS resummation formalism can be adapted for describing nuclear collisions using the 
leading-twist approach discussed in the next section. Since the hard scattering contributions 
are identical in proton and nuclear collisions in the leading-twist approximation, it suffices to 
replace the proton PDFs f a /h{ x ^F) in all parts of the resummed cross section (Q]) by nuclear 
PDFs f a /A( x ^F) for each initial-state nucleus. The smearing function W NP (b, Q) may also 
depend on the type of the nuclear projectile, for example, due to nuclear broadening effects [30] . 
and can be modified by adjusting the a parameter. 

Following this approach, the resummed predictions can be obtained for any set of nuclear 
PDFs that is available [SJ El [TU1 E] ■ In the next Section, we describe the FGS nuclear PDFs 
that will be used in our numerical simulations. 

3 Leading-twist nuclear shadowing in nuclear parton distribu- 
tions 

3.1 Theory roundup 

Nuclear shadowing in hadron-nucleus and real photon-nucleus scattering is a firmly established 
experimental phenomenon. At high energies, the scattering cross section on a nuclear target is 
smaller than the sum of the cross sections on the individual nucleons comprising the nuclear 
target. In the nucleus rest frame, nuclear shadowing is explained as being an example of quan- 
tum mechanical destructive interference between the scattering amplitudes corresponding to the 
interaction with one, two, three, etc., nucleons [31]. The modern theory of nuclear shadowing 
is based on the connection between nuclear shadowing and diffraction [32J , which can be under- 
stood [33] as a manifestation of unitarity reflected in the Abramovsky-Gribov-Kancheli (AGK) 
cutting rules [34j . The accuracy of the resulting theory of nuclear shadowing is rather high, with 
corrections at the level of no more than a few percent. 

Turning to hard processes with nuclei, it was observed that the connection between shadowing 
and diffraction is also valid for any processes initiated by a hard probe and, in particular, for 
deep inelastic scattering (DIS) with nuclei. Combining this connection with QCD factorization 
theorems for inclusive and diffractive DIS, nuclear parton distribution functions (PDFs) for 
small x (x < 0.05) can be calculated. The approach based on this connection is called the 
leading-twist theory of nuclear shadowing [331 E3 [36] and its predictions for nuclear PDFs 
will be used in the numerical analysis in this paper. 
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Figure 1: Multiple scattering series for nuclear quark PDFs. Graphs a, b, and c correspond 
to the interaction with one, two, and three nucleons, respectively. Graph a gives the impulse 
approximation; graphs b and c contribute to the shadowing correction; the interaction with more 
than three nucleons is not shown but included in the final answer. 



The modifications of nuclear parton distributions due to nuclear shadowing are given by a 
multiple scattering series, where each term corresponds to the interaction with one, two, three, 
etc. nucleons. In the graphical form, an example of such a series for the quark distributions in 
nuclei is presented in Fig. [TJ The corresponding expression for the nuclear PDF of flavor j at a 
certain initial scale Qq reads |36t 15]: 

xfj/A&Qo) = Axf j/N {x,Ql) 

f~\ ' \2 /*0.1 p poo poo 

- MA{A-l)Ke \~ lT1 l / ' dxjp/3/f (4) (/3,Q2 )X7P)tmin ) / d % / dzi / dz2 

J- T" T] J x J J— oo J z\ 

x PA (b,z 1 )p A (b,z 2 )e i ^ z ^ m -e-^ il ' lv ^^ , (6) 

where ff^ is the diffractive parton distribution of the nucleon; p A is the nuclear matter 
density; rj is the ratio of the real to imaginary parts of the elementary diffractive amplitude, 
r] = ?R.eA dlS /9m J 4 dlff ~ 0.17. The diffractive PDF depends on two light-cone fractions 

xjp = (M^ + Q 2 )/(W 2 + Q 2 ) and /3 = x/xp, and the invariant momentum transfer t, where 
W is the invariant virtual photon- nucleon energy, W 2 = {q+p) 2 , and M\ is the invariant mass 
squared of the produced intermediate diffractive state, denoted as "X" in Fig.[TJ The longitudinal 
(collinear with the direction of the photon momentum) coordinates z\ and z 2 and the transverse 
coordinate (impact parameter) b refer to the two interacting nucleons; m n is the nucleon mass. 
The t dependence of ff^ can be safely neglected as compared to the strong fall-off of the 
nuclear form-factor for A > 4. As a result, ff enters Eq. ([6]) at t m ™ ~ — x 2 m^(l + M^/Q 2 ) 2 , 
and all nucleons enter with the same impact parameter b. 
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In Fig.[T]and in Eq. ([6]), modification of nuclear PDFs due to the interaction with two nucleons 
is calculated in a model- independent way through the nucleon diffr active PDFs, similarly to the 
case of hadron-deuteron scattering in the Gribov- Glauber theory of shadowing in soft processes. 
The contribution to nuclear shadowing from the interactions with n > 3 nucleons requires 
additional model-dependent considerations, since the interaction of a hard probe (virtual photon) 
with n > 3 nucleons is sensitive to fine details of the diffractive dynamics. Noticing that the 
analysis of diffraction in DIS implies that hadronic fluctuations of the virtual photon have 
predominantly large sizes, one can reliably parameterize the strength of the interaction with 
n > 3 nucleons by a single effective hadron-like cross section of ft . It was argued in Ref. [5] 
that the magnitude of of ft can be reliably estimated using either the color dipole model or 
the model for hadronic fluctuations of the pion. As a result, using these two approaches one 
effectively obtains the lower and upper boundaries on of ft , which in turn correspond to the 
upper (FGS1CLH) and lower (FGS1CLL) boundaries on the predicted nuclear shadowing. 

Note that Eq. © defines xfj/^(x, Qq) only for the shadowing region of x < 0.1 and at Qo = 
2 GeV [5]. In the 0.03 < x < 0.2 interval, there was introduced an additional effect of nuclear 
antishadowing (enhancement) in the gluon channel and modeled by requiring conservation of 
the momentum sum rule. Finally, while it was assumed that in the gluon and sea quark channels 
x fj/A(, x iQo) = Axfj/N(%,Qo) for x > 0.2, the valence quarks are subject to all four types of 
nuclear modifications briefly discussed in the introduction. As a consequence of the DGLAP Q 
evolution from Qq to Q > Qo, the large- x nuclear effects, notably, the EMC effect, start also 
affecting the gluon and sea quark nuclear PDFs, see Figs. 2 and 3. 

3.2 Practical implementation 

We write a nuclear PDF as 

f j/A {xiQ 2 ) = Rj(x,Q 2 ) ■ (Af j/N (x,Q 2 )) , (7) 

where Rj(x,Q 2 ) is the ratio given by the leading- twist theory of nuclear corrections; fjw is the 
free nucleon PDF. For j = u,u,d, and d we have 

Af j/N (x, Q 2 ) = Zf j/p (x, Q 2 ) + {A- Z)f j/n (x, Q 2 ) , (8) 

where A is the atomic number, Z is the number of protons, A — Z is the number of neutrons. 
The subscripts p and n indicate a proton and a neutron, respectively. The PDFs for the neutron 
for u and d (anti-)quarks are obtained via charge symmetry, f u / n (x,Q 2 ) = fd/ P (x,Q 2 ), etc. For 
heavier quarks and gluons, we have 

fj/N(x,Q 2 ) = fj/ p {x,Q 2 ) , j = s,s,c,c,b,b,g. (9) 
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Examples of resulting predictions of the leading-twist theory of nuclear corrections are pre- 
sented in Figs. [2] and [3] , where the ratios of the next-to-leading order (NLO) nuclear ( 208 Pb) to 
nucleon PDFs, fj/A{ x i Q 2 )/[Afj/N( x > Q 2 )]-, ls plotted as functions of Bjorken x at Q = 8 and 90 
GeV, the representative virtualities of the 7* jZ boson encountered in our numerical calculation. 
The two sets of curves (labeled FGS1CLH and FGS1CLL) correspond to the two scenarios for a 3 ft 
(nuclear shadowing) that we have mentioned above. At Q values shown in the figures, one pre- 
dicts significant medium modifications of nuclear PDFs for each parton flavor: large shadowing 
(up to 30-40%) for x < 10~ 2 , antishadowing at 0.05 < x < 0.15, and the suppression originating 
from the EMC effect for x > 0.3 values. The spread between the solid and dotted curves in 
Figs. [2] and [3] is the theoretical uncertainty of the present approach to nuclear PDFs; it is less 
than 10% for A ~ 200 and is much smaller for light nuclei. 

3.3 Experimental tests in the Drell-Yan process 

Next, we wish to outline how the x dependence of the nuclear correction can be investigated 
experimentally by taking advantage of the basic kinematics of the Drell-Yan process. The 
hadronic cross section a for Drell-Yan pair production involves convolutions over longitudinal 
momentum fractions £1,^2 of the hard scattering cross section a with the PDFs, 

da ^ f 1 f 1 da 



dQ 2 dy dQ 



2 



f f da 

f 1 #1 / d6^i,6)4(--lV--l 



^2 



(10) 



"i,2 'Ms.l- v2> " 1 

In the second line, the Dirac delta function is separated from all other terms, indicated collec- 
tively by /i(£i, £2)- It arises due to energy- momentum conservation and constrains the integration 
variables £1 and £2 to lie on a hyperbolic integration contour satisfying K\K 2 = Q^/M^, where 
Ki, 2 = £1,2/2:1,2 - 1, si,2 = (M T /VS)exp(±y), and M T = ^Q 2 + Q 2 T . 

Although £1 and £2 are integrated over an extended range, a large part of the cross section is 
contributed from around the point K\ = K2 = Qt/Mt, where both £1 and £2 are small enough, 
and the corresponding PDFs are large. At this point, 

£2 « x 2c = r e ffe~ y , (11) 



where 



_Mt + Qt , . 

Teff = ~ 7 =— (12) 
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Figure 2: Predictions of the leading- twist theory of nuclear corrections for 
fj/A( x i Q)/[^4/j/Af( x ) Q 2 )] f° r different partons in g2 8 Pb, with Q = 8 GeV, and for the 
indicated parton flavors. The two sets of curves correspond to the two scenarios of nuclear 
shadowing (see the text). 
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Figure 3: Same as Fig. El for Q = 90 GeV. 
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Far away from this point, one of the PDFs is suppressed as £ — > 1, and the cross section is 
reduced. 
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Figure 4: Representative momentum fractions x\ c and X2 C in initial-state particles 1 and 2 as 
functions of y and T e //. 



The pair {xi c ,X2 C } can thus characterize typical values of the momentum fractions in the 
initial-state hadrons 1 and 2 in most of contributing events for the given boson's 4-momentum. 
They can be evaluated from r e ff and y (or, equivalently, y/S, Q, Qt, and y) with the help of 
Eq. (fTTI) . A contour plot for conversion of {r e //,y} into {xi c ,X2 C } is presented in Fig. HI Here, 
the contours indicate constant values of T e ff and y; the corresponding x\ c and X2 C values can be 
read off the horizontal and vertical axis. 

For the following discussion, we also include Table [U listing the ranges of r e ff for the com- 
binations of \/S, Q and Qt values typical for our figures. 

We will argue in Sec. [4] that the behavior of the nuclear correction to the Z/*y* cross section 
can be understood qualitatively from the plots of the nuclear correction to the PDFs in Figs. [2] 
and El taking the x\ c and xi c as the reference momentum fractions. The relation is accurate 
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Table 1: Ranges of r e ff values for kinematical regions shown in the figures. 

within 5-10% and works better when small £ values dominate, since ^(^1,^2) hi Eq. (110p is a 
monotonous function for small £12- More detailed predictions require to know the behavior of 
MCii^) along the full integration range in Eq. (|10p . For the bulk of Z bosons, Qt is much 
smaller than Q, and Eq. (|lip reduces to the familiar 

6,2 ~ xi,2 = (Q/VS)<&. (13) 

However, Eq. (jllj) also clarifies the non-trivial Qt dependence of our numerical results, to which 
we turn next. 

4 Numerical predictions 
4.1 Setup 

The FGS PDFs were implemented in the CSS resummation codes Legacy and ResBos [161117] 
and used to compute Drell-Yan pair production cross sections in two characteristic mass ranges 
of Q = 70 — 110 and 5 — 20 GeV and for the boson rapidity \y\ < 3.5 corresponding to x > 10~ 5 
at all considered energies. ResBos performs resummation of the initial-state QCD radiation 
and decay of massive bosons along the lines described in Sec. [2J At Q « Mz, we computed 
the dominant process of resonant Z boson contribution, as well as small contributions from the 
continuous 7* background and Z — 7* interference. At Q = 5 — 20 GeV, we included the virtual 
photon contribution only, as the Z-mediated cross section is negligible in this mass range. 

Nuclear PDFs are generated according to Eq. ([7]) using a CT10 NLO family of PDFs of the 
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proton [37] • To evaluate the nuclear correction to the differential cross section, we introduce the 
ratio 

/ dajA+B-^ZX) x 

/ \ \ Hv ) with nucl.corr. a\ 

r ^ V ) = t da{A + B^ZX) , • ( 14 ) 

v dv ' no nucl.corr. 

where v = y or Qt, in select ranges of Q and/or y. The denominator is calculated assuming the 
absence of nuclear modifications of nuclear PDFs, i.e., according to Eq. (J7J with Rj(x, Q 2 ) = 1. 
The cross sections are computed for the following values of the center-of-mass energy at the 
LHC: \f S p a = 4.4 TeV and 8.8 TeV in the proton-lead scattering and ^/Saa = 2.76 TeV and 
5.52 TeV in the lead-lead scattering. 

In the perturbative resummed form factor ([2]), the functions A, B, and Cj/ a were evaluated to 
orders , a 2 ., and a s , respectively |144 [38j [39l HQ] . The scale parameters in W pert were chosen 
as {Ci = 6*3 = 26o,C2 = 2}. The nonperturbative contributions were introduced according 
to Ref. [28] using the 6* convention with b max = 1.5 GeV -1 and W NP (b) = exp [— a6 2 ]. We 
take a = 1.1 GeV 2 by default in Z production and a = 0.3 or 1.1 GeV 2 in low-Q Drell-Yan 
production. We assume the same a value in the numerator and denominator of the ratio r a , so 
that the dependence on a cancels well inside r a . The predictions for r a that will be shown are 
hardly sensitive to the value of a that is assumed. 

The Y term was estimated to NLO in photon-mediated subprocesses, and to NNLO in the 
pure Z cross section by using the two- loop correction from Ref. |41| . The renormalization and 
factorization scales were set to 2Q in the Y piece. All these settings are in good agreement with 
ATLAS Qt distributions for pp — > Z/~f*X [32], as has been found in Ref. [28J. 

4.2 Proton-lead collisions 

The r a ratios for proton-lead collisions are shown in Figs. [5H7J for Z boson production in a 
representative interval 70 < Q < 110 GeV and in Figs. [BlfTOl for low-Q Drell-Yan process at 
5 < Q < 20 GeV. The dashed purple line and short-dashed black line were computed using the 
FGS10_H and FGS10_L nuclear PDFs, respectively. 

4.2.1 Z pole region 

We start with the plots of the r a ratio for the Z rapidity distribution in Fig. as they are the 
simplest. In pA collisions, the nuclear correction acts on the PDF of the lead nucleus but not 
on the proton PDF. The nuclear correction depends on which is of order X2 C = i~ e ff exp(-y) 
for the most part. We can therefore understand the overall magnitude of the nuclear correction 
within 5-10% by estimating x<i c for the given y and Qt, as has been argued in Sec. 13.31 
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Since Qt is much smaller than Q for the majority of Z events, a scan of r a (y) over y in 
Fig. [5] essentially translates into a scan over the momentum fraction x in the PDF nuclear 
correction in Fig. [3) Negative (positive) values of y correspond to large (small) values of x = 
(Q/ V S p a) exp(— y). The typical x value (equated to X2 C ) can be found for each y from Fig. d] 
and Table [JJ 

In Fig. Owe see the small-x shadowing at x < 2 • 10 -2 , antishadowing at 2 • 10~ 2 < x < 0.3, 
and large- x suppression at x > 0.3. The exact boundaries between the three regions depend on 
the parton flavor. In the upper inset of Fig. [5l these regions translate into r a < 1 at y > 0, 
r a > 1 at —2.6 < y < 0, and r a < 1 for ?/ < —2.6. The suppression is the most pronounced (of 
order 15-20%) in the forward rapidity regions. The FGS10_L and FGS10_H predictions differ in 
the small-x/large-y shadowing region, but are very close otherwise. 

In the lower inset (obtained for wSpA = 8.8 TeV), we observe a qualitatively similar behavior 
of r a as in the upper inset (for S p a = 4.4 TeV), but the antishadowing "hill" is shifted toward 
lower y values, —3.2 < y < —0.8, as a result of reduced r e // for the same Q. 

The measurement of Qt of the lepton pair provides an additional handle for probing the x 
dependence. For y fixed, a larger Qt value leads to a larger r e ff and X2 C as compared to Qt = 0. 
Ranges of r e // in the Qt intervals shown in the figures are listed in Table [TJ 

We also observe that the flavor dependence of scattering contributions is distinct at small 
and large Qt- At Qt *C Q, the scattering proceeds largely through qq annihilation, while at 
Qt w Q the qg Compton scattering becomes competitive. The analysis of Qt dependence can 
thus identify nuclear corrections attributed to the gluon PDF. 

The Qt dependent ratios r CT (Qy) for pA — > Z production are plotted at 4.4 and 8.8 TeV 
in Figs. [6] and [7] in the specified bins of y. The Qt bins were chosen as in the pp —> ZX 
measurements by ATLAS [42]. The rapidity dependence is not symmetric with respect to y = 
in the pA case. Hence the Qt distributions are plotted separately for positive and negative y. 

In both figures, the value of r a at low Qt in each bin is in rough correspondence with 
the r a (y) value for the corresponding y in the rapidity distribution. This is expected, since 
most Z events that contribute to the rapidity distribution have small Qt- As Qt increases, the 
r a (QT) changes non-trivially, especiall in the positive y bins (right columns in the figures), where 
significant small-x shadowing (of order 15-20%) at Qt — > is superceded by antishadowing of 
up to 10% at Q T > 50 - 100 GeV. 
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4.2.2 Low-Q region 

In the low invariant mass range, 5 < Q < 20 GeV, Z-mediated e + e~ production is negligible 
compared to virtual photon contributions. In this case Table [TJ and Fig. U]tell us that the Drell- 
Yan process probes much smaller values of X2, C; of order a few 10 -3 . The r a rapidity distribution 
in Fig. [5] shows that small- x nuclear shadowing dominates across the full rapidity range, with the 
exception of a small band of large negative rapidities below —2.5. The magnitude of shadowing 
exhibits a monotonous increase from zero at y ~ —2.5 to 25-30% at y ~ +3. 

Even when shadowing dominates across in some y interval, antishadowing may still occur at 
this y if Qt is high enough. This is illustrated by Figs. [9] and [101 showing Qt dependence of 
the r a ratio in the low-mass range. In the bins with y < in the left columns of Figs. Oland fTUl 
r a is smaller than 1 (larger than 1) at Qt = {Qt = 20 — 30 GeV). For y < —3, even large- x 
suppression of up to 10% occurs at the largest Qt values close to 100 GeV. 

At positive y in Figs. [91 and [TUl (right columns), small Qt shadowing is more prominent. 
In this y region, corresponding to x < 2 • 10 -3 for Qt = 0, shadowing can reach 40% in 
the lowest Qt bins. As the transverse momentum increases above Qt = 10 GeV, r a remains 
below unity, in consistency with the x^ c value that is still quite small. Detailed variations of r a 
are now more intricate, as a wide range of the momentum fractions £2 contributes at these y: 
£,2,min < £2 < 1, where £,2,min ~ ^2 is much smaller than 0.01. As a result, r a at Qt > 10 GeV 
shows both the dependence on the small- x shadowing model (from £2 < 0.01) and a "hill" from 
the antishadowing and enhancements in the valence PDFs at £2 > 0.01. 

It can be further shown that the qg Compton scattering contributes of order 80% at the 
largest Qt shown in the figures. The Drell-Yan process at low Q and large Qt is therefore an 
excellent probe of nuclear corrections to the nuclear PDFs, in full analogy with pp case |43l I44| . 

4.3 Numerical results for lead-lead collisions 
4.3.1 Z pole region 

The nuclear correction is more complex in lead-lead collisions, since both incoming beams con- 
tribute to the nuclear effects. In this case the rapidity distribution is symmetric, as can be seen 
in Fig. [TTJ The r e // values with Qt = are equal approximately to 0.033 at ^/Saa = 2.76 
TeV and 0.016 at \/Saa = 5.52 TeV and correspond to mild antishadowing and small-x shad- 
owing, respectively. Consequently the y = region in Fig. [11] exhibits antishadowing of 2-3% 
for y/ Saa = 2.76 TeV (upper inset) and 5% shadowing for ^/Saa = 5.52 TeV (lower inset). In 
the central region, predictions based on two nuclear models are almost indistinguishable. 
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As \y\ increases, the strong small- x shadowing in one of the initial-state nuclei overtakes 
the mild increase in r a due to antishadowing associated with the second nuclei. At the largest 
rapidity shown (|y| — > 3.5) the shadowing reaches up to 30% and has different magnitude in the 
two nuclear models. 

Qt dependence of r a is illustrated in Figs. [T2l and PT3l in three bins of \y\. We observe the 
strong small- x shadowing (15-20%) at the smallest Qt values, which is replaced by antishadow- 
ing (up to 10%) at intermidiate Qt, and eventually by large-x suppression at the highest Qt 
attainable. The shadowing is generally stronger for \/Saa = 5.52 TeV (Fig. [13]) as a result of 
smaller typical r e ff, of about 0.001 for Qt = and 0.036 for Qt = 100 GeV. 

4.3.2 Low-Q region 

In the low invariant mass regime illustrated by Figs. [THI161 the shadowing is most pronounced; 
antishadowing in one nuclei never overcomes shadowing in the counterpart nuclei in either y 
distribution or Qt distribution. Antishadowing produces characteristic bumplike features in the 
lower inset of the y distribution in Fig. Q3] and Qt distributions in Figs. [15] and [T6] but never 
overcomes the overall suppression that results in r a < 1. In the Qt distribution, the shadowing 
is generally stronger at the smallest Qt (up to 40%) than in the y distribution (25-30%). The 
profile of r a at the intermediate and large Qt value again displays a "hill" , similarly to the pA 
case. 

4.4 Nonperturbative smearing of the transverse momentum distribution 

One of the key uncertainties in determination of the nuclear PDF corrections from Drell-Yan 
Qt distributions is associated with the small-Qr nonperturbative function W NP (b,Q) that is 
not completely known in nuclear scattering. In our analysis, W NP (b, Q) has been determined 
from pp collisions (cf. Sec. [2]) and parametrized as 

W NP (b, Q) = exp {-a(Q)b 2 ) , (15) 

with a(Q) = 1.1 GeV 2 at Q « Mz and a = 0.3 — 1.1 GeV 2 at low Q. In the previous subsections, 
the ratio r a has quantified the effect of the nuclear correction to the nuclear PDF and was 
constructed to be almost independent of a. But for the actual Qt distributions, the nuclear 
modifications may also arise from the nonperturbative function WnpQ>,Q), which has impact 
on production of Drell-Yan pairs at Qt < 20 GeV. 

Analysis of the experimental data on Drell-Yan dilepton production in pA collisions at = 
40 GeV finds [45J a very small ^-dependent broadening of the Qt distribution. At this energy 
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(corresponding to x > 0.1) and Qt < 3 GeV, A{Q T ) grows approximately linearly with A 1 / 3 
and reaches 0.1 GeV 2 for the heaviest nuclei. Moreover the analysis indicates that the data at 
lower energies are consistent with the energy-independent A(Q T ). Physically, the broadening of 
the Qt distribution can be interpreted as multiple rescattering of the incoming quark off A 1 / 3 
nucleons. At smaller momentum fractions £, the broadening from rescattering may be larger 
due to the increase of the gluon density at small x. At the same time the A dependence of the 
broadening effect may become weaker due to the leading-twist nuclear shadowing. 

In the CSS resummation formalism, broadening is accommodated by modifying the nonper- 
turbative kr distribution in the nucleus at small x in the region of nuclear shadowing. The 
nuclear shadowing is stronger for small virtualities, so that effective k T at scale 1 GeV may be 
somewhat larger than in the nucleon case. 

To account for this possibility, and the fact that a has uncertainties even in pp collisions, 
we give the plots of da /cIQt for a rather wide range of a in a representative case of lead- lead 
collisions at \/Saa = 2.76 TeV in the shadowing model FGS10_H. 

In Fig. [T7] we compare predictions for da/dQT at 5 < Q < 20 GeV for the range a = 
0.3, 0.7, 1.1., 1.5, and 1.9 GeV 2 , divided by da/dQ T for a = 0.3 GeV 2 . The value of a = 0.3 GeV 2 
is of order of the average a value in pp collisions in this Q interval in fixed-target (large-x) Drcll- 
Yan production [18] . Larger a values would reflect the increased nonperturbative kr in the nuclei 
and result in broader Qt distributions that may be more typical for small x values [46J and/or 
nuclear targets. In this specific example, increasing a from 0.3 to 1.9 GeV 2 increases the average 
(Q T ) in the affected interval < Q T < 20 GeV by 2.5 GeV 2 (from 28.6 to 31.1 GeV 2 , which 
is much larger than A(Q T ) in a fixed-target experiment discussed above). The dependence of 
{Qt) on a reflects the interplay of perturbative and nonperturbative QCD contributions and 
varies with y/S and Q. 

Similarly, in Fig. [18] we show predictions for da/dQT at 70 < Q < 110 GeV for the range 
a = 0.7, 1.1, 1.5, and 1.9 GeV 2 , divided by da/dQT for the nominal a = 1.1 GeV 2 . The range 
0.7 — 1.5 GeV 2 is of order of the current 95% experimental uncertainty in a(Q) in pp — > ZX 
|27] . The curve for 1.9 GeV 2 is an example of more extreme broadening than in pp collisions. 
Note that Qt distributions in Z production are generally less sensitive to the nonperturbative 
smearing than in low-Q Drell-Yan process. 

In the two figures, we observe that the strongest dependence on a(Q) occurs at Qt < 4 GeV, 
where the cross section varies between 80 and 150% (85 and 130%) in the low-Q (high-Q) region. 
The variations are much milder (±10% and ±5%) at 4 < Qt < 20 GeV and practically vanish 
at Qt > 20 GeV. We see that, while the nonperturbative smearing is relevant at Qt below 5 
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GeV, it leaves the Qt distributions intact at larger Qt and hence will not influence sensitivity 
to the nuclear PDF corrections in most of the Qt range. 

5 Conclusions 

We studied production of the neutral Z and 7* gauge bosons with their subsequent decay into 
lepton pairs in proton-lead and lead-collisions at the LHC and examined the role of nuclear 
medium modifications of nuclear PDFs. By extending the ResBos code for the Collins-Soper- 
Sterman Qt resummation to the case of nuclear parton distributions, we analyzed the transverse 
momentum Qt and rapidity y dependence of the production cross section. We examined r a , 
the ratio of the differential production cross sections with and without nuclear modifications of 
quark and gluon PDFs in nuclei, and found unambiguous correspondence between the predicted 
behavior of r a and the pattern of nuclear modifications of nuclear PDFs. At central or moderate 
rapidities, there may be a region with r a > 1 as a consequence of antishad owing of the valence 
quark distributions in nuclei. For forward and backward rapidities, we generally predict that 
r a < 1 due to the suppression of the quark and gluon distributions in nuclei by nuclear shadowing. 
For large rapidities the ratio r a can also discriminate between different scenarios (magnitudes) 
of nuclear corrections, while the large-Qy region probes the gluon nuclear correction. A variety 
of measurements that can be carried out with pA and AA Drell-Yan production make it an 
informative test of the QCD evolution in high-energy heavy-nuclei scattering and of the pattern 
of medium modifications of nuclear PDFs. 

The ResBos code and resummed grids for simulations in this study can be downloaded from 
http:/ /hep. pa.msu.edu/resum/." 
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Figure 5: The nuclear correction ratio r a (y), as denned in Eq. fH]) . plotted for 70 < Q < 110 
GeV vs. the lepton pair rapidity y in proton- lead collisions. 
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Figure 6: The nuclear correction ratio r a (QT) for the Qt distribution and 70 < Q < 110 GeV 
in proton-lead collisions at ^S p a = 4.4 TeV, in six bins of y. 
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Figure 7: Same as Fig. EJ for ^~S^ = 8.8 TeV. 
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Figure 8: The nuclear correction ratio r a (y), as denned in Eq. (fTl"|) . plotted for 5 < Q < 20 GeV 
vs. the lepton pair rapidity y in proton-lead collisions. 
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Figure 9: The nuclear correction ratio r a (QT) for the Qt distribution and 70 < Q < 110 GeV 
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Figure 10: Same as Fig. H for ^/S^ = 8.8 TeV. 
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Figure 12: The nuclear correction ratio r a (QT) for the Qt distribution and 70 < Q < 110 GeV 
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Figure 13: Same as Fig. [12j for \/SpbPb 
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Figure 14: The nuclear correction ratio r a (y), as defined in Eq. (fl4"|) . plotted for 5 < Q < 20 
GeV vs. y in lead-lead collisions 
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Figure 15: The nuclear correction ratio r CT ((5r) for the Qt distribution and 5 < Q < 20 GeV in 
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Figure 16: Same as Fig. QSJ for ^5^1 = 5.52 TeV. 
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Figure 17: Dependence on nonperturbative smearing of the Qt distribution in lead-lead collisions 
for 5 < Q < 20 GeV at s/S^A = 2.76 TeV. 
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Figure 18: Same as Fig. El for 70 < Q < 110 GeV. 
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